Nernst Effect as a Signature of Quantum Fluctuations in Quasi-ID Superconductors 
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We study a model for the transverse thermoelectric response due to quantum superconducting 
fluctuations in a two-leg Josephson ladder, subject to a perpendicular magnetic field B and a 
transverse temperature gradient. The off-diagonal Peltier coefficient (a xy ) and the Nernst effect are 
evaluated as functions of B and the temperature T. The Nernst effect is found to exhibit a prominent 
peak close to the superconductor-insulator transition (SIT), which becomes progressively enhanced 
at low T. In addition, we derive a relation to diamagnetic response: a xy = —M/To, where M is the 
equilibrium magnetization and To a plasma energy in the superconducting legs. 
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In low-dimensional superconducting (SC) systems 
(thin films, wires and Josephson arrays), fluctuations 
of the SC order parameter field lead to broadening of 
the transition to the SC state, and give rise to anoma- 
lous transport properties in the adjacent normal phase 
PP. While close to or above the critical temperature T c 
thermally excited fluctuations dominate these conduction 
anomalies, quantum fluctuations are expected to domi- 
nate at low temperatures T <C T c , where superconductiv- 
ity is weakened due to, e.g., the effect of a magnetic field, 
disorder or repulsive Coulomb interactions. Their most 
dramatic manifestation is the onset of a superconductor- 
insulator transition (SIT) when an external parameter 
such as magnetic field or thickness is tuned beyond a 
critical point [H [3] . 

A striking signature of the fluctuations regime, which 
attracted much attention in the recent years, is the 
anomalous enhancement of transverse thermoelectric ef- 
fects in the presence of a perpendicular magnetic field 
B. In particular, a substantial Nernst effect measured 
far above T c , e.g., in the underdoped regime of high-T c 
superconductors [U [5] and disordered thin films [B] . As 
the Nernst signal (a voltage developing in response to 
a temperature gradient in the perpendicular direction) 
is typically small in ordinary metals, its magnification 
in such systems has been attributed to the dynamics of 
thermally excited Gaussian SC fluctuations [7HS], or mo- 
bile vortices above a Kosterlitz-Thouless [TU] transition 
[IT] . Theoretical studies have also been extended to the 
quantum critical regime of SC fluctuations [12]. 

Conceptually, the above mentioned theoretical models 
share a common intuitive idea: in the phase-disordered, 
vortex liquid state, vortex flow generated parallel to a 
thermal gradient (V a T) naturally induces an electric field 
(E x ) in the perpendicular direction. Consequently, the 
general expression for the Nernst coefficient 

^ Ex Pxx^xy Pxy&yy 

(WyT)B ~ B U 

is overwhelmingly dominated by the first term, dictated 
by the off-diagonal Peltier coefficient a xy . The latter is 
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FIG. 1: (color online) A scheme for measurement of the 
Nernst effect in a Josephson ladder subject to a magnetic 
field B perpendicular to the plane, and a temperature differ- 
ence between the top (Ti) and bottom (T2) SC wires. Dashed 
lines represent Josephson coupling. 



an interesting quantity: a transport coefficient, yet inti- 
mately related to thermodynamic properties. In partic- 
ular, it was found to be proportional to the diamagnetic 
response 13 [IT]. In the clean limit, it was shown to 
encode the entropy per carrier [13] [14]. Note, however, 
that the implied contribution to v [Eq. ([!])] is not deter- 
mined by a xy alone, but rather its product with the elec- 
tric resistivity p xx - Observation of a large Nernst signal 
therefore necessitates a reasonably resistive normal state. 

In this paper, we present a theory for the transverse 
thermoelectric coefficients and their relation to diamag- 
netism in the quasi one-dimensional (ID) superconduct- 
ing device depicted in Fig. [Tj in which the geometry 
dictates an appreciable Nernst effect in the fluctuations- 
dominated regime. The device considered is a two-leg 
Josephson ladder subject to a perpendicular magnetic 
field B, where a small temperature difference between the 
legs induces voltage along the ladder due to transport of 
vortices across the junction. At low T, one expects vortex 
transport to be dominated by quantum tunneling. This 
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system serves as a minimal setup for observing trans- 
verse thermoelectric effects; the relative simplicity of the 
model describing the quantum dynamics of SC fluctua- 
tions allows an explicit evaluation of a xy , v and the mag- 
netization density M in a wide range of parameters. In 
particular, we investigate their behavior when the wires 
parameters are tuned through a SIT, and find a promi- 
nent peak in v close to the transition, which becomes 
progressively enhanced at low T. We further confirm the 
proportionality relation between a xy and —M, however 
the pref actor is I /To, with To the plasma energy scale, 
rather than 1/T as found in the 2D case [51 [TTj. 

As indicated in Fig. [I] the system consists of two SC 
wires of length L parallel to the x direction separated by 
a thin insulator layer of width W, which allows a weak 
Josephson coupling J per unit length. In each of the 
separate wires (n = 1,2), the ID quantum dynamics of 
fluctuations in the phase of the SC condensate is governed 
by a Hamiltonian of the form (in units where K = 1) 



v 

2?r 



dx 
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here (j) n {x) is the collective phase field, and 8 n (x) its 
conjugate field (obeying \<j> n {x),d x O n {x')\ = iit5{x' — x)) 
which denotes Cooper pair number fluctuations. This 
model can be viewed as describing, e.g., the contin- 
uum limit of a Josephson chain [3] , where the Josephson 
coupling (Ej) and charging energy (E c ) between adja- 
cent SC grains are related to the parameters of H n by 
K = \J E c /Ej and v = \/E c EjTra, with a the grain size. 
The Josephson coupling between the wires is given by 



Hj = — J / dx cos{</>i — 4>2 — Qx} 



(3) 



where q is the deviation of the vortex density in the junc- 
tion area from a commensurate value: 



n {WB , 1 

= 2n mod - 

V $ a 



(4) 



in which (for K = c = 1) $o = ?r/e is the flux quantum. 

Assuming the hierarchy of scales J a -C T -C v/a (with 
T an average temperature of the system), Hj [Eq. ^] 
can be treated perturbatively. Note that the first inequal- 
ity justifies this approximation for an arbitrary value of 
the Luttinger parameter K in Eq. for K < 2 and 



sufficiently small q, the Josephson term becomes rele- 
vant [15], and induces a SC phase where fluctuations in 
the relative phase between the wires are gapped [121 [17] 
in the T — > limit. In turn, additional perturbations 
such as inter- wire charging energy [161 117j and disorder 
[T5] generate a transition to an insulating T — > phase 
for sufficiently large K. Since, as shown below, in both 
extreme phases the Nernst effect is strongly suppressed, 
we focus our attention on the intermediate T regime. In 
addition to Hj, we introduce a weak backscattering term 
due to disorder of the form 



H D = ^2 dx( n (x) cos{29 n (x)} 



(5) 



71=1,2 



( n (x) = , Cn{x)Cn'(x') = D6(x - x')6 n n < 



where overline denotes disorder averaging. This term 
generates the leading contribution to the resistivity p xx , 
and thus to v via Eq. ([lj. 

We now consider a temperature difference AT = T\ — 
T 2 between the top and bottom wires (see Fig. [T]), each 
assumed to be at equilibrium with a separate heat reser- 
voir, and focus first on the transverse Peltier coefficient 
a xy . Evaluating the electric current I x induced along the 
ladder for AT < T yields a xy = I X /AT. Alternatively, 
Oixv = Ta xy can be derived from the heat current I x h ^ 
generated by a voltage difference V y between the wires. 
We show below that the result of both calculations is the 
same, as dictated by the Onsager relation [8j. 

The electric current is given by the expectation value 



I x = 2eir ( a 



(6) 



where the current operators 6 n (x,t) — i[H,8 n ] (H 
II: + Hj where Hq = H i + H^) are evaluated pertur- 
batively in Hj using the interaction representation. The 
leading contribution to I x arises from the second order: 

9 n (x,t) = U{t)9 ( °\x,t)U\t) where 

eW(x,t) = ~d x <f> n (x,t), (7) 

U(t) = 1 + i$dtxHj{ti) - / dh / dt 2 Hj(t 1 )H J {t 2 ) . 

— oo — oo — OO 

Employing Eq. (|3| and inserting the resulting expres- 
sions for 9 n in Eq. mh, we obtain 



Trev J 



J dt\ J dt 2 J dx\ J dx 2 sin [q(xx — x 2 )] i^m 



e 2 f Ti 



Ftj (xi— x 2 ;ti— *2)g— it Ft 2 (xi-x 2 ;t 1 -t 2 ) 



[{d x F Tl (x - xx;t - ti) + d x F Tl {x 2 -x;t 2 - t)} - {d x F T2 (x -x\;t- h) + d x F T2 (x 2 -x;t 2 - t)} 
- {d x F Tl (x - x 2 ;t - h) + dxF Tl (xi -x;t 2 - t)}+ {d x F T2 (x - x 2 ;t - ti) + d x F T2 {xi -x;t 2 - t)}] 



(8) 
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where we use the Boson correlation function Fr(x;t) = ([</> n (x , i) — (fi n (0, 0)] 2 ) at fixed temperature T [15] : 



F r (a:;t) = -log 



2 

-^57= r (sinh [ttT ( - - i + i-sign(i) ) sinh \ttT \-+t- i-sign(i) ) X 
'K z a i ± z I L \u v )\ I \v v /J J 



For a small temperature difference between the wires (AT «T = ^p-) and L — > oo, this yields 
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(10) 



where B(z, w) is the Beta function. 

We next consider the alternative setup where an elec- 
tric voltage V y is imposed between the top and bottom 
wires (at uniform T), and evaluate the heat current I. 
induced along the ladder. For J = 0, the local heat cur- 
rent operator is given by [19] 



(h) 



(11) 



n=l,2 



The voltage bias corresponds to a difference in chemical 
potentials in the two legs, /ii,2 = ±eV y , which introduce 
constant shifts of d x B\,2 by ±ireV y /vK. The heat current 



Ah) 



[/(i)j£ 0) [/t(t)\ [with U(t) expanded to second 



order in Hj as in Eq. fffc], is hence given by 
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FIG. 2: (color online) Isotherms of v as a function of K for 
v/27va — IK, T = 0.06, 0.07, 0.08, 0.09, 0.1K. 



j(>0 



envVy 
K 



(U(t)(d x (f>x ~ d x <f> 2 )W(i)) . (12) where the constant prefactor v cx LDJ 2 and 



The resulting expression coincides with (V y T/AT)I x , 
with I x given by Eq. j8l) for AT < T. We thus con- 



firm that a xy = I. 



(h) , 



To derive the Nernst coefficient, we employ Eq. ([!]) 
noting that within our level of approximations, a yy and 
p xy vanish due to particle-hole symmetry. The Nernst 
signal in the setup depicted in Fig. [T] defined as v = 
V/ATB, is hence determined by the product of a xy [Eq. 
(10)] and the longitudinal resistance of the ladder R xx . 
To leading order in H D [Eq. ^] [3j Q22 [18], 
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(13) 

At low magnetic field such that q — 2irWB/$ <C T/v, 
this yields an expression for v sa a xy R xx / B of the form 



vaT{K) 



(2-koT 



K+4-6 



(14) 



_r 2 (f)r(i-f)r(i){^(f)-V/(i-f)} 



2 2 W(l-f)F(f)F(i 



(15) 

[r(z), ip 1 ' [z) are the Gamma and Trigamma functions, re- 
spectively]. The dependence of v on K, the parameter 
which tunes the SIT in the SC wires, is depicted in Fig. 
[5] for different values of T <C v/a. In this regime, v ex- 
hibits a pronounced maximum at K*(T), slightly below 
the transition from SC to insulator (K c = 2 [31120]). As 
T is lowered, the peak becomes progressively enhanced 
and K* ~ y2 as dictated by the rightmost exponential 
factor in Eq. (14 1. This non- monotonous behavior can 



be traced back to the competition between electric resis- 
tance (which signifies the rate of phase-slips), and a xy 
(which signifies the strength of diamagnetic response). 

We next investigate the relation of a xy to the diamag- 
netic magnetization density M = —(l/LW)(dF/dB), as- 
suming that the ^-dependence of the free energy F is 
restricted to the flux in the junction area, i.e. the pa- 
rameter q in Eq. pi). To leading order in Hj and at low 
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magnetic field B = (<E> /27rJU)<2, 

2 T 

Sj — —J J dx J dr cos {<fii{x, r) 



The resulting expression for M is identical to Eq. ( 10 ) 



up to a constant prefactor. We thus obtain the relation 



M 
To 



To — — — V E c Ej 
ira 



(17) 



where the last equality associates the energy scale To with 
the plasma energy in the Josephson chain forming the 
legs of the ladder. This resembles the linkage pointed 
out in earlier literature, except the thermal energy scale 
T in the prefactor is replaced here by the characteristic 
scale of quantum dynamical phase-slips. 

We finally note that the remarkable relation to the 
entropy per carrier a xy ~ —(s/B), derived for clean sys- 
tems [T31 [H], does not hold here. Indeed, this relation 
can be recovered, e.g., employing the Boltzmann equa- 
tion for an ordinary conductor in the limit lo c t — > oo 
(with u c the cyclotron frequency and r a mean free time 
[2"T]). In contrast, for uj c t <C 1 the same calculation yields 
OLxy ~ Bt 2 . In our case, the latter limit is appropriate: 
while translational invariance holds in the x-direction, 
charge conductance along the y-direction is governed by 
weak tunneling between two discrete points, ~ J. We 
hence expect a xy ~ BJ 2 , in accord with Eq. (10 1. 



To summarize, we studied the transverse thermoelec- 
tric coefficients due to quantum SC fluctuations in a 
Josephson two-leg ladder, and their relation to diamag- 
netism. Most importantly, we predict a large Nernst sig- 
nal, particularly at low T where a pronounced peak is 
predicted close to the SIT. This behavior reflects a sub- 
tle interplay between diamagnetism (favored in the SC 
phase), and dynamical phase- fluctuations (which prolif- 
erate in the insulator). As a concluding remark, we ex- 
pect a qualitatively similar effect to hold in 2D SC films 
(or an infinite stack of such ladders) . Possibly, it can also 
explain some properties of the existing data [6 . 
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